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ABSTRACT 



The light front analysis of 7r~ mesons in He(Li,C), C-Ne, C-Cu and O-Pb collisions 
is carried out. The phase space of secondary pions is divided into two parts in one of 
which the thermal equilibrium assumption seems to be in a good agreement with the 
data. Corresponding temperatures T are extracted and their dependence on (Ap * A^) 1 / 2 
is studied. The results are compared with the predictions of the Quark-Gluon String 
Model (QGSM). The QGSM satisfactorily reproduces the experimental data for light and 
intermediate-mass nuclei. 
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equilibrium, characterized by the temperature T. Analysis in light-front variables. Com- 
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1. INTRODUCTION 



In the present paper we continue the study of tt~ mesons produced in the relativistic 
nucleus-nucleus collisions in terms of the light front variables. The choice of the light front 
variables is due to the fact, that these variables seem to be more sensitive to the dynamics 
of interaction as compared to the well-known Feynman variables xf and rapidity Y . 

In the previous publication [1] we performed the light front analysis of n~ mesons 
produced in Mg-Mg collisions at 4.3 GeV/c. On the basis of this analysis we were able to 
separate the region of the phase space, where the thermal equilibrium seems to be reached. 
It is interesting in this connection to perform the same analysis for other pairs of nuclei 
( He(Li,C), C-Ne, C-Cu, O-Pb) in order to study the dependence of the corresponding 
characteristics on the atomic number of colliding pairs of nucleus . The data were obtained 
on the SKM-200-GIBS facility of the Dubna Joint Institute for Nuclear Research. 

2. EXPERIMENT 

SKM-200-GIBS consists of a 2m streamer chamber, placed in a magnetic field of ~ 0.8 
T and a triggering system. The streamer chamber was exposed by beams of He , C , O , 
Ne and Mg nuclei accelerated in the synchrophasotron up to a momentum of 4.5 GeV/c 
per incident nucleon. The solid targets in the form of thin discs with thickness 0.2^-0.4 
g/cm 2 ( for Li the thickness was 1.59 g/cm 2 and for Mg 1.56 g/ cm 2 ) were mounted within 
the fiducial volume of the chamber. Neon gas filling of the chamber served also as a 
nuclear target. The triggering system allowed the selection of "inelastic" and "central" 
collisions. 

The central trigger was selecting events with no charged and neutral projectile specta- 
tor fragments (with P/Z > 3 GeV/c ) within a cone of half angle Q c h, Q n = 2° or 3°. The 
trigger mode for each exposure is defined as T (Q ch ,6 n ). Thus T(0,0) corresponds to all 
inelastic collisions. For inelastic collisions He-Li and He-C all charged secondaries were 
measured and central subsamples T(2,0) were selected. The details of the experimental 
setup and the logic of the triggering systems are presented in [1,2]. 

Primary results of scanning and measurements were biased due to several experimental 
effects and appropriate corrections were introduced. The biases and correction procedures 
were discussed in detail in [2,3]. Average measurement errors of the momentum and 
production angle determination for tt~ mesons were < AP/P >= 5%, A© =0.5° for 
He-Li, He-C, C-Ne, C-Cu, O-Pb. 
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3. THE DATA ANALYSIS IN TERMS OF THE LIGHT FRONT 

VARIABLES 



The light front variables ^ in the centre of mass frame for the inclusive reaction 
a + b — > c + X are defined as follows [1]: 



t± t E±p z ,E+\p. 



where s is the usual Mandelstam variable, E = yp 2 + p\ + m 2 and p z are the energy and 
the z - component of the momentum of produced particle. The z -axis is taken to be 
the collision axis i.e. p z = p%. The upper sign in Eq. (1) is used for the right hand side 
hemisphere and the lower sign for the left hand side hemisphere. In order to enlarge the 
scale in the region of small it is convenient also to introduce the variables 

C ± = Tlnl^l (2) 

The invariant differential cross section in terms of these variables looks as follows: 

E do_ = da = 1 da 

dp 7r d^dpj, it d^dp^ 

The light front variables have been introduced by Dirac [4] and they are widely used 
now in the treatment of many theoretical problems (see, e.g. original and review papers 
[5-10]). 

The analysis has been carried out for the 7i~ mesons from He(Li,C), C-Ne, C-Cu and 
O-Pb collisions in the nucleus-nucleus centre of mass system. 

The number of events for all pairs of nuclei and corresponding trigger modes are 
listed in Table 1. Due to the small statistics and average multiplicities, the data of He- 
Li and He-C collisions has been united and thus He(Li,C) represents this sample of the 
data. In Fig. 1 the ^ - distributions of 7r~ mesons from He(Li,C), C-Ne, C-Cu and 
O-Pb interactions are presented. These distributions are similar for all analysed pairs of 
nuclei. The principal differences of ^ distributions as compared to the corresponding ij? 
- distributions (Fig. 2) are the following: (1) existence of some forbidden region around 
the point ^ = 0; (2) existence of maxima at some in the region of relatively small 
1^1- 

The experimental data for invariant distributions • dN/dC^ are shown in Fig. 3. 

The curves are the result of the polynomial approximation of the experimental distribu- 
tions. The maxima at are also observed in the invariant distributions • dN/dC^. 
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However, the region > goes over to the region < IC^ and vice versa (see 
Eqs. (1) and (2)). The value of maxima are observed at C ± = 2.0 ± 0.1 for all pairs of 
nuclei. The ^ is the function of the energy (see Eqs. (1), (2)) and does not depend on 
the mass numbers of the projectile (A P ) and target (A T ). 

In order to study the nature of these maxima we have divided the phase space into 
two regions > |C ± | (C += 2.0) and < |C ± | an d studied the p\ and the angular 
distributions of 7r~ mesons in these regions separetely. The numbers of pions in these 
two regions are approximately equal. For example in C-Cu interactions in the region 
IC^l > IC^l the number of pions is equal to -1987 and in < 1(^1 — 2212. In Figs. 4, 5 
the °p\ and the angular distributions of ir~ mesons from He(Li,C), C-Ne, C-Cu and O-Pb 
interactions in different regions of ( + ( ( + > ( + and ( + < ( + ) in the forward hemisphere 
are presented. 

One can see from Figs. 4, 5, that the p\ and the angular distributions of 7r" mesons 
differ significantly in ( + > ( + and ( + < ( + regions. The angular distribution of pions in 
the region ( + < ( + (Figs. 5.b and 5.c) is sharply anisotropic in contrast to the almost flat 
distribution in the region ( + > (Figs. 5. a and 5.c). The flat behaviour of the angular 
distribution allows one to think that one observes a partial thermal equilibrium [11] in 
the region > (1^1 < l^l) of phase space. The slopes of j)\ - distributions differ 
greatly in different regions of C ± (Fig. 4). Thus the values of C± are the boundaries of the 
two regions with significantly different characteristics of n~ mesons. The validity of this 
statement can be seen from the momentum distributions of n~ mesons in the Laboratory 
frame. Fig. 6 presents the momentum disribution of pions from C-Cu collisions in the 
laboratory frame. The shaded area corresponds to the region of ( + > ( + and the non- 
shaded one to the region of ( + < ( + . One can see from the Fig. 6, that these two 
regions almost do not overlap in the momentum space unlike to the cms case (overlap 
~ 45%). The pions from the region ( + > ( + have small momentum, approximately up 
to 0.6 GeV/c as compared to the pions from ( + < (+ (the momentum of pions ranges 
from ~ 0.6 GeV/c to 3 GeV/c). The similar results have been also obtained for the other 
pairs of nuclei. Fig. 7 presents the dependence of < P >i a b on Qi a b for all analysed pairs 
of nuclei (He(Li,C) and C-Ne data are presented with the same symbol because of the 
similarity of their dependences) and Mg-Mg interactions in the ( + > ( + and ( + < ( + 
regions. The shapes of these dependences are different in two regions of ( + . The curves 
are the result of polynomial approximation. < P >i ab decreases and < O >i a t, increases 
with the increasing of Ap, At- 

To describe the spectra in the region ( + > ( + the Boltzmann 

f(E) ~ e~ E / T 
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distribution have been used. 

The distributions 1/n ■ dN/d( + , dN/dp^, dN/dcosQ look in this region as follows 

1 dN fP 2 T , max 



7T d( + 

dN rp 



z . rn ax 



dp\ Jo 
dN fP 



d cos 9 Jo 



/ Ef(E)dp 2 T (4) 
Jo 

f(E)dp z (5) 

f(E)p 2 dp (6) 



max 



where: 



E = ^p 2 + ml, p 2 =pl+p 2 T (7) 

P\max = {( + ^s) 2 - ml 
p z , max = \p\ + m 2 - ((+^S) 2 ]/(-2(+^7s) 
Pmax = {-£ + ^cosQ + yf (£ + t/s) 2 - mlsin 2 Q) I sin 2 Q 

The experimental distributions in the region ( + > ( + have been fitted by the expres- 
sions (4), (5), (6), respectively. The results of the joint fit of the distributions l/n-dN/d( + , 
dN/dp^, dN/dcosQ are given in Table 1 and Figs. 4, 5, 8. They show a rather good agree- 
ment with experiment. In Table 1 the values of the parameter T obtained by fitting the 
data with Boltzmann distribution are presented. In order to determine how the charac- 
teristics vary the analysis have been carried out also for £+=1.9 and 2.1. The results are 
similar, but the joint fit of the distributions is better for (+=2.0 (presented on Figures). 

The spectra of n~ mesons in the region ( + > ( + are satisfactorily described by the 
formulae which follow from the thermal equilibrium. The same formulae when extrapo- 
lated to the region ( + < ( + (Fig. 8) deviate significantly from the data. Therefore in the 
region ( + < ( + the p\ - distributions has been fitted by the formula 



dN 



dp T 

and the ( + - distributions by the formula 

1 dN 



2 ~ a ■ e~^ p T + (!-«)• e-te* (8) 



71 d(- 



l-£+) n = (l- e -IC + l)« (9) 



which is an analogue of the (1 — xp) n dependence - the result of the well-known quark- 
parton model consideration (see, e.g. [12]), which for n~ mesons gives the value n=3. The 
dependence (1 — e _ ^ + l)™ is in good agreement with experiment in the region ( + < (+ anc [ 
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deviates from it in the region ( + > ( + ( Fig. 8). The results of the fit are given in Table 
1 and Figs. 4 and 8. 

Thus in the C ± (^ ± ) distributions we have singled out points C^ 1 (^ ± ) which seperate in 
the phase space two groups of particles with significantly different characteristics. There 
are no such points in the xf and Y - distributions. 

In this paper the Quark Gluon String Model (QGSM) [13] is used for a comparison 
with experimental data. The QGSM is based on the Regge and string phenomenology 
of particle production in inelastic binary hadron collisions [14]. The QGSM simplifies 
the nuclear effects (neglects the potential interactions between hadrons, Pauli blocking, 
coalescence of nucleons and etc.) and concentrates on hadron rescatterings. The QGSM 
includes only low lying vector mesons and baryons with spin 3/2, mostly A (3/2, 3/2) 
resonances. A detailed description and comparison of the QGSM with experimental data 
in a wide energy range can be found in Ref. [1,15,16]. 

We have generated He(Li,C), C-Ne, C-Cu, O-Pb interactions using Monte-Carlo gen- 
erator COLLI, based on the QGSM. The events had been traced through the detector 
and trigger filter. 

In the generator COLLI there are two possibilities to generate events: 1) at not fixed 
impact parameter b and 2) at fixed b. The events have been generated for b. From the 
impact parameter distribution we obtained the mean value of < b >. For the obtained 
value of < b >, we have generated a total samples of A P — A T events. The numbers of 
generated events for all analysed pairs of nuclei are listed in Table I. The two regimes 
are consistent and it seems, that in our experiment the following values of b are most 
probable: b=1.55 fm for He(Li,C); b=2.20 fm for C-Ne; b=2.75 fm for C-Cu; b=3.75 fm 
for O-Pb. 

The experimental results have been compared with the predictions of the QGSM for 
the above mentioned values of b and satisfactory agreement between the experimental 
data and the model have been found. In Figs l.b and 3.b the ^ and - distributions 
of 7r~ mesons from the QGSM calculations are presented together with the experimental 
ones for C-Cu interactions. One can see, that the QGSM reproduces these distributions 
well. The similar results have been obtained for all analysed pairs of nuclei. The QGSM 
also reproduces the p\ and cos© distributions (Figs. 4.c and 5.c). The QGSM data 
show the similar characteristics in the different regions of ( as experimental ones: sharply 
anisotropic angular distributions in the region ( + < ( + and the almost flat distribution 
in the region ( + > ( + ; the slopes of p\ - distributions differ greatly in different regions of 
( + ; the momentum distributions of pions in the laboratory frame in different regions of 
( + have the similar different shape of spectra as experimental ones (Fig 6). Momentum 
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distributions of QGSM data reproduces the corresponding experimental spactra in both 
regions of ( + . The distributions obtained by the QGSM in the region ( + > ( + have been 
fitted by the expressions (4), (5), (6). The results of the fit are given in Table 1 and Figs. 
4.c, 5.c. In the region ( + < ( + the p\ and the ( + - distributions have been fitted by the 
formulae (8) and (9), respectively. The results of the fit are given in Table 1 and Fig. 4.c. 
One can see from the Table 1, that the values of the T extracted from the experimental 
and QGSM data coincide within the errors. The QGSM does not reproduce satisfactorily 
the O-Pb data. This is may be caused by the fact, that QGSM simplifies the nuclear 
effects, which are more pronounced for heavy nuclei. In Ref. [17] it has been indicated 
that the model can be improved by including higher mass baryon resonances and taking 
into account a possible increase of the pion absorption cross section, <jaa^nn in dense 
baryon medium, in comparison with the cross section, obtained from detailed balance 
relation. 

In Fig. 9 the dependence of the parameter T from the Table 1 on (A P * At) 1 / 2 , 
obtained from the experimental and QGSM data, is presented. The temperature for Mg- 
Mg interactions is extracted from [1]. One can see, that T decreases linearly with the 
increasing of (Ap* Ar) 1//2 i.e with the increasing number of participating nucleons. Similar 
behaviour is predicted by the QGSM. 

In our previous article [18] the temperatures of pions in He-Li, He-C, C-Ne, C-Cu 
and O-Pb interactions were obtained be means of inclusive kinetic energy and transverse 
momentum spectra in central rapidity interval (0.5 - 2.1 for light nuclei and 0.1 - 1.8 
for heavy ones), which corresponds to the pionization region and with the c.m.s. angles 
90° ±10°. The pion spectra for He-Li, He-C and C-Ne have been fitted by one exponent and 
for C-Cu and O-Pb by a sum of two exponents, or two temperatures Ti and T 2 (describing 
the low and high momentum part of the spectrum). The temperatures extracted by the 
light front analysis for light pairs of nuclei are less about (15-20) % as compared to those 
obtained in Ref. [18]. For heavy pairs of nuclei the temperatures are more close to the low 
temperature T\. It seems obvious, as the thermal equilibrium region corresponds to lower 
momenta. It should be mentioned that the extraction procedures of T in the light-front 
variables and in Ref. [18] are quite different and it seems, that different regions of phase 
space are seperated by these two methods. 

It is interesting to compare the temperatures extracted by means of light front analysis 
with those obtained in the GSI experiments (FOPI, KAON and TAPS- Collaborations, 
see, e.g. [19,20]). The T in the GSI experiments have been obtained in a same manner as 
in our Ref. [18]. The numerical values of the parameter T for pions in Au-Au collisions 
at 1 A GeV and our values for heaviest colliding pair are close to each other. It seems 
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interesting in this connection to perform the light front analysis of the GSI, AGS and SPS 
data. 



4. CONCLUSIONS 



The light front analysis of n~ - mesons in the relativistic nucleus-nucleus collisions is 
carried out. The results of this paper confirm the conclusions of our previous publication 
[1], that the phase space of secondary n~ mesons is divided into two parts, in one of 
which the thermal equilibrium seems to be reached. Corresponding temperatues are 
obtained from the fitting of the data by the Boltzmann distribution. The characteristics 
of the 7r~ mesons (the momentum, angular, p\ - distributions) in these two regions differ 
significantly. The dependence of T on (A P * A T ) 1 ^ 2 has been studied. The temperature 
decreases with increase of (A P * A?) 1 ^ 2 '■ The experimental results have been compared 
with the QGSM. The model seems to be in a good agreement with the data excluding 
O-Pb collisions - the heaviest colliding pair. 
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FIGURE CAPTIONS 



Fig. 1 The ^ - distribution of tt~ mesons a) from * - He(Li,C), A - C-Ne, o - 
O-Pb interactions, b) From C-Cu interactions: o - the experimental data, A - the QGSM 
data. 

The curves are the result of polynomial approximation of the experimental data. 

Fig. 2 The x F - distribution of tt~ mesons from * - He(Li,C), A - C-Cu, o - O-Pb 
interactions. 

Fig. 3 The C ± - distribution of n~ mesons a) from * - He(Li,C), A - C-Ne, o 
- O-Pb interactions, b) From C-Cu interactions: o - the experimental data, A - the 
QGSM data. 

The curves are the result of polynomial approximation of the experimental data. 

Fig. 4 The p\ distribution of 7r~ mesons from * - He(Li,C), A - C-Ne, o - O-Pb 
interactions a) for ( + > ( + . b) For ( + < £+. c) From C-Cu interactions for ( + > ( + 
o - experimental and A - the QGSM data; for ( + < □ - experimental and * - the 
QGSM data. 

The solid lines - fit of the experimental data in the regions ( + > and ( + < ( + by the 
Eqs.(5) and (8), correspondingly. 

Fig. 5 The cos© distribution of n~ mesons from * - He(Li,C), A - C-Ne, o - O-Pb 
interactions a) for ( + > ( + . b) For ( + < £+. c) From C-Cu interactions for ( + > (+ 
o - experimental and A - the QGSM data; for ( + < ( + , □ - experimental and * - the 
QGSM data. 

The solid lines - fit of the experimental data in the region ( + > ( + by the Eq. (6) and in 
the ( + < ( + by the polynom. 

Fig. 6 The momentum distribution of mesons from C-Cu interactions in the 
laboratory system. The shaded area corresponds to the region of ( + > ( + . 

Fig. 7 The dependence of < P >i ab on Q iab in the regions ( + > ( + (bottom data) 
and C + < C + ( to P data ) for □ - He(Li,C) and C-Ne, o - C-Cu, * - Mg-Mg, o - O-Pb. 
The curves - result of polynomial approximation. 
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Fig. 8 The (l/V) • dN/d( + distribution of ir~ mesons from C-Cu interactions, o - 
experimental data, the solid line - fit of the experimental data in the region ( + > ( + by 
the Eq.(4), the dashed line - fit of the experimental data in the region ( + < ( + by the 
Eq.(9). 

Fig. 9 The dependence of the parameter T on (Ap * A T ) l l 2 for He(Li,C), C-Ne, Mg- 
Mg, C-Cu and O-Pb. o - the experimental data, A - the QGSM data. The dashed line 
is a result of linear approximation of the experimental data. 



TABLE CAPTIONS 



Table 1. Number of the events, trigger and the results of the joint fit of the 
distributions 1/n ■ dN/d( + , dN/dp^, dN/dcosQ of n~ - mesons by Eqs. (4), (5), (6) in 
the region ( + > ( + and l/ir ■ dN/d( + distributions by Eq.(9) in the region ( + < ( + . 
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Table 1. Number of the events, trigger and the results of the joint fit of the distribu- 
tions 1/V • dN/d( + , dN/dp T , dN/dcosQ of 7r~ - mesons by Eq. (4), (5), (6) in the region 
( + > ( + and 1/ir • dN/d( + distributions by Eq.(9) in the region ( + < ( + . 
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14 



120 




160 




Figure 1: The ^ - distribution of 7r _ mesons a) from * - He(Li,C), A - C-Ne, o - O-Pb 
interactions, b) From C-Cu interactions: o - the experimental data, A - the QGSM data. 
The curves are the result of polynomial approximation of the experimental data. 
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Figure 2: The Xf - distribution of n mesons from * 
He(Li,C), A - C-Cu, o - O-Pb interactions. 
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Figure 3: The ^ - distribution of n~ mesons a) from * - He(Li,C), A - C-Ne, o - O-Pb 
interactions, b) From C-Cu interactions: o - the experimental data, A - the QGSM data. 
The curves are the result of polynomial approximation of the experimental data. 
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Figure 4: The distribution of 7r~ mesons from * - He(Li,C), A - C-Ne, o - O-Pb 
interactions a) for ( + > ( + . b) For ( + < ( + . c) From C-Cu interactions for ( + > ( + 
o - experimental and A - the QGSM data; for ( + < ( + □ - experimental and * - the 
QGSM data. The solid lines - fit of the experimental data in the regions ( + > ( + and 
( + < ( + by the Eqs.(5) and (8), correspondingly. 
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Figure 5: The cos® distribution of 7r~ mesons from * - He(Li,C), A - C-Ne, o - O-Pb 
interactions a) for ( + > ( + . b) For ( + < ( + . c) From C-Cu interactions for ( + > ( + 
o - experimental and A - the QGSM data; for ( + < ( + □ - experimental and * - the 
QGSM data. The solid lines - fit of the experimental data in the region ( + > ( + by the 
Eq.(6) and in the ( + < ( + by the polynom. 
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Figure 6: The momentum distribution of tt mesons from C-Cu interactions in the labo- 
ratory system. The shaded area corresponds to the region of ( + > ( + . 
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Figure 7: The dependence of < P >i a b on iab in the regions ( + > 
C+ (bottom data) and ( + < ( + (top data) for □ - He(Li,C) and 
C-Ne, o - C-Cu, * - Mg-Mg, o - O-Pb. The curves - result of 
polynomial approximation. 
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Figure 8: The (1/vr) • dN/d( + distribution of n~ mesons from C- 
Cu interactions, o - experimental data, the solid line - fit of the 
experimental data in the region ( + > by the Eq.(4), the dashed 
line - fit of the experimental data in the region ( + < by the 
Eq.(9). 



22 



100 



He-Li(C) 

C-Ne 



Mg-Mg 



50- 



o 



C-Cu 



> 

CD 



O-Pb 



60- 



40 



O (~\ I I I I I I 

zu 10 20 30 40 50 60 

(A P *A T ) 1/2 

Figure 9: The dependence of the parameter T on (A P * At) 1 / 2 for 
He(Li,C), C-Ne, Mg-Mg, C-Cu and O-Pb. o - the experimental 
data, A - the QGSM data. The dashed line is a result of linear 
approximation. 
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